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Abstract 

In this paper a convergent series expansion is constructed to solve the pre- 
scribed mean curvature equation ±V •(Vu/y 7 ! ± |Vii| 2 ) = nH for n-dimensional 
hypersurfaces in R n+1 (+ sign) and R 1 '™ (— sign) which are graphs {(x,u(x)) : 
x E R n } of a smooth function u : R n — > R, and whose mean curvature function 
H E Cq (R n , R) is not too large in norm and integrable. Our approach is inspired 
by the Maxwell-Born-Infeld theory of electromagnetism in R 1 ' 3 , for which our 
method yields the first systematic way of explicitly computing the electrostatic 
potential u for regular charge densities p oc H and small Born parameter. There- 
fore, after the general n-dimensional problems have been treated with the help of 
nonlinear Hodge theory and Banach algebra estimates, our approach is reworked 
in more detail for n = 3. 

©2010 The copyright for this preprint resides with the authors. Its reproduction, in 
its entirety, for non-commercial purposes is permitted. 
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1 Introduction 



In this paper we develop a constructive approach to the prescribed mean-curvature 
problem for hypersurfaces which are graphs of some scalar function u over the entire 
W 1 . While we use only standard tools from the theory of elliptic PDE, we believe 
that we combine these in a novel efficient way to find solutions to this so-called non- 
parametric prescribed mean-curvature problem. The scalar second-order problem for 
u is reformulated in terms of an equivalent first-order Hodge problem. It is then 
shown that for suitably small mean curvature functions the Hodge problem admits a 
convergent power series expansions which can be solved term-by-term through explicit 
quadratures. The convergence of this series in a Banach algebra can be proven with 
the help of complex analysis. 

Our approach grew out of our attempts to solve the nonlinear Maxwell-Born-Infeld 
equations of the electromagnetic fields in R 1 ' 3 , i.e. when n = 3. In the electrostatic 
limit, recalled in Appendix A, the electrostatic potential is proportional to the function 
u, and the charge density is proportional to the mean-curvature function. Thus our 
method also provides a systematic way of finding non-special electrostatic solutions 
to the Maxwell-Born-Infeld field equations when the charge density is suitably small. 
However, the smallness condition for H misses the important case of electrostatic 
fields generated by point charges [TJ EJ [2QJ [2H E5j [23], [13] , corresponding to maximal 
hypersurfaces with conical singularities [IB] . 

In the next section we begin our quest by stating the scalar prescribed mean- 
curvature problem for graphs over R n in both Euclidean and Minkowskian geometry. 
In section 3, we write down a formal power series solution to the scalar problem(s) for 
integrable prescribed curvature functions which are small in a Cq norm; the rest of our 
paper is concerned with studying the convergence of the formal power series solution in 
all dimensions n G N. In section 4 we use nonlinear Hodge theory to first reformulate 
the scalar second-order PDE problem into a first-order system for two 1-forms, then 
state the corresponding formal power series solution for these nonlinear Hodge systems 
when the prescribed curvature function is small, and finally prove its convergence 
using Banach algebra arguments and a small amount of elementary complex analysis. 
Because of its special relevance to the mathematical theory of electromagnetism we also 
reformulate the three-dimensional problems in terms of equivalent vector problems, 
and also present some additional details; this is done in section 5. In section 6 we 
indicate some generalizations of our approach to related quasi-linear PDE problems in 
divergence form. 



2 



2 The scalar prescribed mean-curvature equation 



The mean curvature function H of a hypersurface in R n+1 (n = 1, 2, 3, ...) which is the 
graph {(x,u(x)) : x G R n } of a real-valued C 2 function u can be computed from u as 

H = n-^-^==, (1) 



u 



(see, e.g., [30], and the appendix to section 14 in [21]), while the analogous result for 
space-like hypersurfaces in Minkowski space R ,n reads 

(see, e.g., [JJ El [221 12E] ) ■ The inverse problem is to find a function u, satisfying certain 
asymptotic conditions, whose graph over R n describes a hypersurface in R n+1 , respec- 
tively R 1 '™, for which a prescribed function H : R n — > R is its mean curvature. In that 
case ([!]) and ([2]) become second-order quasi-linear elliptic PDEs for the unknown u, 
which we lump together as 

Vn 

±V ■ . = nH, (3) 

'1± |Vd 2 V ; 



the upper sign for the Euclidean, the lower for the Minkowskian setting. 

Regarding the problem of hypersurfaces which are graphs over R n , we note that 
the scalar prescribed mean-curvature equation can be solved in great generality when 
n — 1 (by simple quadratures, obviously). And for n = 2 powerful complex variable 
techniques can be put to work, see [TBI 1221 12S1 El EH [23 E], in particular when 
H = 0. However, for n > 3 no general solution method exists. While one finds many 
special results and general (non-)existence and regularity results, there is not much in 
terms of explicitly constructing non-special solutions — except when H is small. For 
small H one can in principle compute solutions by running the iterations of the fixed 
point maps of the Schauder theory, or use the familiar variational methods; see, e.g. 
[3] for a recent contribution. 

In the remaining sections we develop an alternative, non-variational, more direct 
constructive PDE approach which produces entire graphs in all dimensions nGN for 
quite general mean curvature functions H as long as they are small. Our approach 
treats the Euclidean and Minkowskian problems on an equal footing. 
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While we here are interested in entire solutions over M. n motivated by the electro- 
magnetic Born-Infeld model, and also by some problems in the theory of spacetime 
structure (e.g. [TJ [3j [TH EH]), with some mild modifications our approach should be 
adaptable to the prescribed mean-curvature equation in bounded domains with small 
Dirichlet data for u (e.g. [301 El EH E5]); this then also includes the minimal / maxi- 
mal hypersurface problem (H = 0) with small Dirichlet data (e.g. [TB I 1^7 ] 128 1 \T7 \ 154]). 
The parametric prescribed mean-curvature problem (e.g. [381 M H2l EH SQl EH]), 
which also captures embedded hypersurfaces which are not graphs over MJ 1 and even 
non-embedded hypersurfaces, has a different structure, however. 



3 Formal solution for small prescribed mean-curvature 

Now let H G Cq n L 1 be small in the following sense: H(x) = eH (x), where H (x) 
is a fixed function, and e G M+ can be chosen as small as we please. Then one can 
expand the scalar prescribed mean-curvature equation fl3]) into a hierarchy of (linear) 
Poisson equations all of which can be solved explicitly in the usual way. The problem 
is to prove the convergence of the so-obtained formal power series (in e) solution. 

Thus, for small but finite e the analyticity properties of the nonlinearity of the 
scalar prescribed mean-curvature equation suggest to seek a solution u in form of the 
power series Ansatz 



u\x) 



X>V p) (*), (4) 



with each independent of e. Inserting this Ansatz into the equation ([3]) and sorting 
according to powers of e, we find a hierarchy of linear Poisson equations. 

In particular, we find recursively that for each even p, i.e. p = 2k for k G N, «( p=2fc ) 
satisfies 

Au^ k) = 0, (5) 

so that for general k G N is an entire harmonic function. We also find that 
satisfies 

±Au {1) = nH , (6) 
which is solved by G Cq'" given by 

u m {x) = T^ I | MV l i Vy + hPix) (7) 
< a / \x — v\ 
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when n^2, and by 

u^(x) = Tj§r l [ Ho^ln^— d'y + h^ix) (8) 

J |*£ y\ 

when n = 2. In ([7]) and ([H]), ftA 1 ^ is an arbitrary entire harmonic function. Carrying on 
we also find recursively that for odd p > 1, i.e. p = 2k + 1 for A; e N, the functions 
u (2fc+i) sa tj s fy linear Poisson PDE 

±An (2fc+l) = y . ^(2^+1) ; (g) 

where [/( 2fc+1 ) is a polynomial in the Vit^ with odd £ < 2fc + 1, viz. 

^(3*+!) = W 2(fc -" )+1) ^Mf ^ W 2 ^- 1+1) -W 2 ^ +1) , (10) 

/i=l j=l W 3 j= i=l 

h-j 

K 

where \£\k = ^2 h-, the £i take any non- negative integer values, and where 
i=i 

Mf = (±iy K J , j " (ii) 

is the j-th Maclaurin coefficient of l/y/l T z (with Mjf := 1). Assuming enough 
regularity of U^ 2k+l \ the general regular solution to the linear Poisson PDE OH]) is 
given by 



when n^2, and by 



u {2k+1 Hx) 



+ fceN (12) 

T^n /"v-t/ (2fc+1) (l/)ln^^d 2 y + /i( 2fc+1 )(x), fc e N (13) 
1 1 J \x — y\ 

when n = 2; here, h^ 2k+1 ^ is an arbitrary entire harmonic function. 

To remove the ambiguity in the choices of harmonic functions h&\ in this paper 
we impose the asymptotic condition on u that 

u m (x) x T ^ J H (y)d n yj-^- 2 (14) 



when n^2, and by 



W «(x)x T r|r| jH (y)d 2 y\n^ (15) 

when n = 2. As a consequence, we find that h^' = for all p EN. 

It follows that u = j2T=o e2k+lu(2k+1) i with uW S iven b Y © and «( 2fc+1 ) by (JI2||. 
recursively for fceN when n^2, respectively given by their two-dimensional analogs 
when n = 2, and with vanishing harmonic functions h^' = for all p E N, is a formal 
series solution for our u problem Q, with H = eH . To prove the convergence of this 
formal power series solution we next resort to nonlinear Hodge theory. 



4 Equivalent Hodge problems 

For all dimensions n E N the second-order prescribed mean-curvature equation for 
either sign is equivalent to a nonlinear Hodge system of first order for some 1-form. 
There are two mutually dual formulations. 

4.1 The first-order Hodge systems 

Consider a 1-form oo satisfying 

du = 0; (16) 



±5^JL== = nH. (17) 
1 ± M 



Then this first-order system is equivalent to the second-order equation (J3]) by identi- 
fying co = Vu ■ dx. 

Dual to the above system is the following. Defining a 1-form 

T = ± vW (18) 

which can be inverted to yield 

" = ± /i~~i — 12 ' ^ 

v 1 =f M 2 
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we see that the above first-order equations are equivalent to 



5r = nH, (20) 



r 



d—=== = 0. (21) 



irl 2 



Clearly, f l2Tj) implies that there is scalar a such that 

T 



IT r 



da. (22) 



Up to a sign and an additive constant, a = u, of course. 

We end this subsection with the observation that we can also carry out the dif- 
ferentiation in (12ip and multiply with (1 =p |r| 2 ) 3 / 2 to find the cubic version of (12T1) . 



= (1 :p \r\ 2 )dT±d \\t\ 2 At. (23) 
We shall come back to this apparently simplest formulation in our next-to-last section. 

4.2 The small curvature r-hierarchy 

For small H = eHo, with Hq e Cq fl L 1 , we can recycle our series solution Ansatz for 
either one of the two dual first-order Hodge problems, and solving one leads also to 
the solution of the other. Yet we found it more convenient to work with the r problem 
in its square root version. Thus we set 

oo 

r = ^eM p) , (24) 
p=i 

with each independent of e. Inserting this Ansatz into the pair of equations (12"0I) . 
( 12TI) . and sorting according to powers of e, we find a hierarchy of linear equations. 
In particular, satisfies 

<5r (1) = nH , (25) 

together with 

dT (1) = 0. (26) 
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Consistency with our asymptotic conditions for u requires that for H E Cg D L 1 we 
select the unique solution E C^'" of the pair of linear equations f )25|) . fl26|) . given 
by the exact 1-form 



r^(x) := -d 'f- 1 ' d"y (27) 



F -2/ 



for n^2; when n = 2 the logarithmic kernel has to be used. Note that for n > 2 we 
even have e Co'". 

Having the exact 1-form (127]) . we next find that r*- 2 -* satisfies 

5r (2) = (28) 

together with 

dr {2) = 0, (29) 

and at infinity. Clearly, the pair of linear equations (J2"8"j) . f[2"9"j) has a unique 

solution r^ 2 ^ which tends — )■ at infinity, given by the trivial solution = 0. 

Carrying on we now find recursively that each t^ p > for even p, i.e. p = 2k for k E N, 
satisfies 

<?r< 2fc > = (30) 
rfr (2fc) = 0, (31) 

with r( 2fc ) at infinity, so that = for general k EN, while for odd p = 2k + 1 
with fc e N, we find the pair of linear first-order PDE 

rfr (2fc+i) = dT (2fc+i) j (33) 

where T( 2fc+1 ) is a polynomial in the with £ < 2k + 1, viz. 

fc ft j 

r (2fc+i) _ _ T (2(fc-fc)+i) J-J r (2f 2l _ 1+ l). r (2<? 2l +l) 5 ( 34 ) 

h=l j=l \e\ 2j =h-j i=l 

where \l\2j and Mj have their previous meaning. The pair of linear first-order PDE 
( )32|) . (|33|) has a unique solution r( 2fc+1 ) E C 6 llQ given by 

r (2fc+i) = pr (2fc+i) j fc G N (35) 
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where P : C£' a — > C b ' a projects onto the co-exact subspace of C\ ,a ' . Since ST^ 2k+1 ^ G 
Cq'° nL 1 , we have explicitly 



for n 7^ 2; when n = 2 the logarithmic kernel has to be used. Note that n > 1, then 
we even have r^ 2fc+1 ' G Cq' q . 

Thus we have formally solved our r problem (|2"U|) . ( 12~T|) with iJ = ei?o in terms of 
the formal series solution r = Y^T=o e 2fc+1 r^ 2fc+1 \ with r^) given by ( 12~T|) and r^ 2k+1 ^ for 
A; G N recursively given by f )36|) . We next investigate its convergence. 

1 Absolute convergence of the formal r series solution 

To show that the formal series r = YlT=o e 2k+l T^ 2k+l ^ converges absolutely, in C b ' a when 
n = 1 and in Cq'" when n > 2, we first have to check that each r^ 2k+1 ^ as given in f )36|) is 
in C^' a , respectively Cq' q . But when n = 1 this follows inductively from the facts that 
T W G C£' a , that C 6 1,a is a Banach algebra, and that P : C b ' a — > C b ' a ; hence, each partial 



sum of r = YlT=o e 2k+1 T^ 2k+1 ^ is in C b ' a . Similarly, when n > 2, then G Cq' q and 
so r^ 2fc+1 - ) G Cq' q and each partial sum of r = YlT=o e 2k+1 T^ 2k+1 ^ is in Cq ,ol . It remains 
to show that the right hand side of the norm estimate ||r|| < Xlfclo e 2fc+1 ||r^ 2fc+1 ^ || 



We now estimate all norms ||r( 2fe+1 )|| for k G N in terms of the 2k + 1-th power of 
Ht^H- For notational convenience we extend the list of T( 2fc+1 )s defined for k G N by 
( |34|) to the case k = by defining 



Now, ||r( 2fe+1 )|| = ||PT( 2fc+1 )|| for k G N, and since P : C l b a -> Cj' a , respectively P : 
Cq° — > Cq°, is a projector, we have the estimate: ||r( 2fc+1 ) || < ||T( 2fc+1 )||. Substituting 
the RHS of fl34|) for T^ 2k+1 \ repeating the standard inequality || J2i T i\\ — 2~2i ll r ill> then 
using the inequality 1 1 Tj- 1 1 < ||rj|| ||tj|| valid in Banach algebras (here C b ' a or Cq' q ), 
and applying repeatedly the identity ||r < - 2a+1 ^|| = ||PT^ 2a+1 )|| followed by the projector 
estimate ||PT( 2a+1 )|| < ||T( 2a+1 )|| for the various pertinent values of a > 1 (no estimate 




(36) 



converges. 



T (1) := r {1) . 



(37) 
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is necessary when a = 0), and using \MJ\ = M- := Mj, for k > 1 we obtain 

k h 2j 

n T (2fe+i)|| < || r (2(fc-fc)+i)|| ^ \mj\ ^2 nii r(2 ^ +i) n 

h=l j=l \i\v=h-j i=l 

k h 2j 

< ^HT^-^II^M, ^ fj||T^ +1 )||. (38) 

/i=l 3=1 K|2j=ft-J »=1 



Setting fc = 1 and recalling the definition ( 157]) we obtain the estimate 

\\T®\\ < §||r«|| 3 (39) 

Now suppose that for all k = 1, A;* there exists some i?2fc+i such that 

|| T (2*+i)|| < J R 2fe+1 || r ( 1 )|| 2fe + 1 . (40) 

Then the estimate ( 158]) guarantees that (140]) is true also for k = fc* + 1, and since 
A;* > 1 is arbitrary in this induction step while (159"]) says that the estimate is true for 
k* = 1, it follows that ( 158]) is true for all G N. Our inductive proof that ( 140]) holds 
for all k G N also yields that i?2fc+i is recursively defined for k G N by 

k h 2j 

-fiWi = -R 2 (fc-fe)+i y jj-r^+i, (41) 

h=l j=l \l\n=h-j i=l 

with i?x := 1. 

Recall that we want to show that ^^L e 2fc + 1 || r ( 2fc+1 ) || < oo for sufficiently small e, 
so since we have shown that ||r( 2fc+1 )|| < ||T( 2fc+1 )|| < i2 2 fc+i||r^|| 2fe+1 , it suffices to show 
that Ylh=o -^2fc+i( e ll T ||) 2fc+1 < oo for sufficiently small Her^H =: ||r 9 ||. Setting now 
||r 9 || =: £, we note that the formal power series G(£) := YlT=o R-2k+it; 2k+1 is nothing 
but the formal generating function of the i?2fc+i, in the usual sense that, formally, 
R2k+i = G^ 2k+1 \0) I (2k + 1)1. So our task is to show that the generating function is 
analytic about £ = with non-zero radius of convergence. 

With the help of the recursion relation (|4"T]) we readily find that is the positive 
inverse function of g t- > £ given by 

Z = 2 9-^== 2 (42) 
V 1 - 9 
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defined for positive £ near £ = 0, with G(0) = 0. Since the function g i->- £ given by 
(j42p is analytic about (7 = (with radius of convergence = 1) and has unit derivative 
at g — 0, there now exists an open neighborhood of £ = on which there is defined 
a unique inverse function £ h- >■ g — which vanishes at £ = 0, has unit derivative 
at £ = 0, and satisfies fl4"2"j) . Thus, in particular, the Maclaurin expansion of 
converges to and it generates the recursion coefficients i^fc+i- 

We now determine the finite radius of convergence of the power series for G(£) 
about £ = 0. Setting G(£) = sin\I/(£) we see that (pf2j) becomes the function ^ i-> £ 
given by 

£ = 2 sin?/; -tan^, (43) 

with £ = when ^ = 0. Since ip H- sin ^ is an entire function, which vanishes for ip = 
and has unit derivative there, the radius of convergence of the Maclaurin series of £ H- 
G(£) = sin \I/(£) coincides with the radius of convergence of the Maclaurin series of £ H- 
\I/(£). This radius of convergence in turn is determined by those £ value(s) closest to £ = 
at which the derivative of ip i— > £ = 2 sin ^ — tan ■?/> vanishes (possibly asymptotically 
should £ — )■ £oo when |^| — >■ oo suitably). But this ip derivative is 2cos , ?/ ; — 1/ cos 2, 0, 
and it vanishes iff 2 cos 3 ip — 1, which gives 2 1//3 cos?/' G {1, e* 2 " 7 ' 3 , e* 4?r//3 }. A calculation 
now gives the radius of convergence of G(£) about £ = as 

& = (2 2 / 3 - 1) 3/2 . (44) 

Thus we have shown that our formal series solution converges absolutely if ||r 9 || < 
( 2 2/3 _ ^ 3 / 2 This com pi e tes our convergence proof. 

We remark that because of our use of the projector estimates we cannot conclude 
that || rj < (2 2 / 3 -l) ' is a necessary criterion for absolute convergence of our formal 
power series solution. Indeed, for any radially symmetric H G Cq the formal power 
series reduces to its first term, all other terms being identically zero, so convergence 
is a trivial issue and holds for any size of e, then. 

4-2.2 Spin-off: series solution of the u and u problems for small H 

The series solution for r immediately yields the series solution for the w-problem as a 
spin-off; namely, 



oj = 

fceN 
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Alternately, one can of course use the r-series solution and then apply (TT9"j) to map r 

to LO. 

Lastly, the series solution for r also yields the series solution for u as a spin-off; 
namely, defining 

oo 

p{x) = J2$T i2k+1 \x) (46) 

k=0 

we find that p 6 Cq D L 1 , and so 

^) = ^/^fpi d ^ ( 47 ) 
for n ^ 2; for n = 2we use the logarithmic kernel. 

5 Prescribed mean-curvature graphs over K 3 

Because of its relevance to the electrostatic Maxwell-Born-Infeld equations, we now 
re-develop our solution theory in slightly more detail for the special dimension n = 3. 
The first step is to recast the prescribed mean-curvature problem for graphs over M. 3 
in vector notation. There are two mutually dual first-order vector versions. 

5.1 The equivalent first-order vector problems 

Let u solve ((3]) for whichever sign and, for the pertinent sign, define the vector field w 
by 

w = ±Vn (48) 
Then, for the pertinent sign, w satisfies 

V--=Z== = 3H (49) 
a/1 ± \w\ z 

together with 

Vxw = 0. (50) 

Given a solution w which vanishes at infinity, a solution u to ([3]) with monopole asymp- 
totics at infinity is obtained by line integration of ( 14H1) . 



u ( x ) = =f / w (y) ■ &y- ( 51 ) 
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The first-order system (14"9I) . (JSUJ) we call the u w problem," the second-order equation 
([3]) is henceforth referred to as the "m problem." 

Dual to this first-order system of vector PDE is an alternate first-order system for 
a vector field v defined for a solution u of ([3]) by 

Wu . , 

(52) 



'I ± |Vm| 2 
The elementary inversion of f|52|) yields 

VM = ± /T^W ( 53 ) 
V 1 =F M 2 

For either sign, v satisfies 

V-v = 3H. (54) 
On the other hand, since Vw is curl-free, taking the curl of ( 1531) gives 

V x = = 0. (55) 



The two vector PDE ( 1541) . ( 1551) also form a closed system of nonlinear first-order PDE 
for v which we call the u v problem." Having solved (154"]) plus ( )55|) for u vanishing at 
infinity, a solution to the prescribed mean-curvature equation(s) ([3]) with monopole 
asymptotics is obtained in terms of a line integral of ( )53|) . 



/■°° ffu) 

u(x) = T . yy > 2 . -dy. (56) 

5.2 Hierarchical vector series solution for small H 

For H e Cq flL 1 small in the sense that if (a;) = eH (x), where H (x) is a fixed function, 
and e G M + as small as we please, our convergent power series solution method for the 
r problem translates into the following power series solution method for the v problem. 
One can in fact work with either one of the two dual first-order vector problems, and 
solving one leads also to the solution of the other, but working with the v problem is 
more convenient than working with the w problem. 
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5.2.1 Convergent series solution of the v problem 

Our power series Ansatz for the nonlinear first-order v problem reads 

oo 

v {x) =^e p v { P\x), (57) 
P =i 

with each v w independent of e, and e small but finite. Inserting this Ansatz into the 
pair of equations ([MD, (JSSD , and sorting according to powers of e, we find a hierarchy 
of linear equations. 

In particular, satisfies 

V ■ v {1) = 3H , (58) 
together with 

V x = 0. (59) 

For H e Cq fl L 1 the pair of linear equations (135]). fl59l) has a unique solution G Cq'" 
given by the gradient field 

v W(x):=-V [ M^MdV (60) 

Having the gradient field (!60|) . we next find that ?/ 2 ) satisfies 

V-t» (2) = (61) 
together with 

V x v {2) = 0. (62) 

Clearly, the pair of linear equations flBTj) . fl6"2"|) is uniquely solved by the trivial solution 
v^=0. 

Carrying on we now find recursively that each v^' for even p, i.e. p = 2k for k E N, 
satisfies 

VV 2fc) = (63) 
V x w (2fc) = 0, (64) 
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so that v ( 2fc ) = for general k £ N, while for odd p > 1, i.e. p = 2fc + 1 for G N, we 
find the pair of linear first-order PDE 

V . v (2fc+i) = o ( 65 ) 
Vx/ + " = Vxl/( a+1 ', (66) 

where V A( - 2fc+1 - ) is a polynomial in the with £ < 2/c + 1, viz. 

fc h j 

y(2k+i) _ _y^ v (2(fc-*)+i) V] n^^i-i+iJ.^^i+i)^ ( 67 ^ 

/i=l i=l l<l2j-= i=i 

where again \£\2j and Mj have their earlier assigned meaning. The pair of linear 
first-order PDE (|65|) . (I66I) has a unique solution given by 

where P: C ' a — >■ Cq'" projects onto the solenoidal subspace of C*o' a ; for V ■ \A 2fc+1 ) e 

c °' a nL\ 

py (2A+1) (x) = y (2*+l) (x ) + V / J*I_: _W d 3 y _ (6Q) 

y F-y| 

So far we have proved that v = YlT=o e 2k+1 v^ 2k+1 \ with v^ 1 ' given by flBTjj) and t>( 2fc+1 ) 
for fceN recursively given by f lBDj) . is a formal series solution for our v problem (j53|) . 
( 155|) with H = eHo. Almost verbatim to the proof of section 3 we can now prove its 
convergence in Cq ,c * for small enough e. Thus the series solution is a classical solution. 

5.2.2 Spin-off: series solution of the w and u problems for small H 

By the duality between the v and w problems it follows right away that as a bonus of the 
solution to the v problem one also obtains the solution to the w problem. Specifically, 
by (|48p and (1531) . if v solves the v problem then the solution to the w problem for the 
pertinent sign reads 

V 1 =F \ v \ 

Inserting the series solution (EH) for v into the RHS of ([70]), recalling that v^ =2k ^ = 
while i>W is given by (16"Uj) and the ti(p =2fc+1 ) by (I68p . (I69p . we find the series 

oo 

w ( x ) = J2^ P w (p) {x), (71) 

p=l 
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with each independent of e and given as follows: 



w M{x) = vM(x), (72) 
and for odd p — 2k + 1, with k G N, 



w^ +1 >(x) = V / ^— r^d 6 y, (73) 

J \x-y\ 

with V^ 2fc+1 )(x) given in terms of the i/ 2fe+1 ) by ( 1671) : for all even p = 2k with k E N, 
we have 

W ( 2n )(x)=0. (74) 

Each term in the w expansion is a gradient field, as it should, by ( 148]) . 

Finally, having the solution to the w problem in form of a series of gradient fields 
converging in C Q ' a , we can invoke (148]) one more time and directly read off the solution 
in Cq'° to the original u problem as the series 

= £ fr^i^ ( 75 ) 

where 

oo 

p(x) = ±^e 2fc+1 V-^ 2fc+1 \ (76) 

fc=0 

with l/( 2fc+1 ^(x) given in terms of the t>( 2fc+1 ) by (IBTj) when fceN, and by (l37j) when 
fc = 0. 

We close this subsection with the remark that the solution to the w problem ob- 
tained here as spin-off of the solution to the v problem does not express u/ 2fc+1 ) as 
an operator applied to some explicit expression involving all lower order w^ k ~ h ^ +1 \ 
h = 1,2, ...,k. We leave it to the interested reader to prove that the series solution 
Ansatz ( !7T|) inserted directly into the PDE of the w problem leads to expressions for the 
w (p) involving only all lower order w^ <p \ which can be converted into the expressions 
involving only all lower order ?/ <p ) given here. 
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5.3 The Helmholtz decomposition 



As shown by Helmholtz, every vector field in R 3 can be decomposed into a sum of a 
divergence-free and a curl-free field. The decomposition is not unique because a har- 
monic field can be transferred from one summand to the other. This non-uniqueness 
can be removed by imposing suitable boundary or asymptotic conditions on the in- 
dividual summands; in our case asymptotic vanishing conditions. The Helmholtz de- 
composition sheds an interesting light on the w and v problems; or rather, on the v 
problem. 

Indeed, for the w problem the Helmholtz decomposition just reverses the steps 
taken by us to arrive at it in the first place. Thus, we write w = w g + w c , where 
w g is a gradient field (hence, curl-free) and w c is a curl of some vector field (hence, 
divergence-free). Clearly, V • w = V • w g and V x w = V x w c , so equation (I50|) 
together with the vanishing asymptotics for w c implies right away that w c = 0, hence 
w = w g , and f H9|) becomes a closed equation for w g = ±Vw. Of course, this is exactly 
the original second-order scalar equation (j3J). 

As to the t>-problem, writing v = v g + v c , where v g is a gradient field (hence, curl- 
free) and v c is a curl of some vector field (hence, divergence- free) , and registering that 
V • v — V • v g and V x v = V x v c , we get 



(together with V • v c = 0). We thereby have obtained a closed set of first-order vector 
equations for the vector field v g separately from v c , plus a conditionally closed set of 
first-order vector equations for the vector field v c , conditioned on v g being given. Both 
sets of PDE are supplemented by the asymptotic conditions that v g and v c vanish at 
spatial infinity; with this asymptotic condition the Helmholtz decomposition becomes 
unique. 

As announced, the Helmholtz decomposition sheds an interesting light on the v 
problem. Since the curl-free component of the Helmholtz decomposition is determined 
completely and independently of the divergence-free Helmholtz component v c of v , 
namely thusly: For H G Cq fl L 1 the linear equation (ITT)) has a unique gradient field 



V ■ v g = 3H 



(77) 



(together with V x v g = 0), and 








(78) 
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solution E C ' a given by 

«,(*):= -V / i^d 3 !/, (79) 

y k — 2/| 

we can reinterpret our power series approach as a method to solve the v c -equation (I78p 
given small v g . Indeed, note that v g = ei;' 1 ' given earlier in (160 j) (note also that with 
n = 3 we have er^ = t> 9 • da;), so we can just separate the v 5 part off from our power 
series solution for v. The remaining power series yields the solution v c to f JTHj) given 
small v g . 

5.4 The cubic nonlinearity of the v problem 

We now come back to our earlier remark that the nonlinearity of the r problem, i.e. here 
the first-order vector problem for v, is effectively cubic. As we will see, this implies some 
computational simplifications (though, interestingly enough, no simplifications for our 
solution theory, yet!). In addition, this "cubic version" of our v problem reveals some 
interesting a priori differential identities which are satisfied by any solution but which 
are obscured by the original formulation of the v problem. 

Carrying out the curl operation in fl55|) and multiplying through with (1 =p \v | 2 ) 3 / 2 
we find that, away from singularities (which occur, for instance, when 1 — \v | 2 — > 0), 
any solution v of ( 155]) satisfies 

= (l =F |f| 2 ) V x v ± |~|V|f | 2 j x v. (80) 

This equation already exhibits a cubic nonlinearity, yet it can be further manipulated 
into a more concise alternate format. 

Namely, using first the identity |V|v | 2 = (v ■ V)f + v x (V x v) well-known from 
vector analysis, we find that flHOl) is equivalent to 

= (l =F M 2 ) V xv±[(v V)v + v x (V x v)) x v. (81) 

We next employ the identity v x (V x v) x v = |w| 2 V x v — v(v- V x v), well-known from 
vector algebra, which for any solution of flHOl) simplifies to v x (V x v) x v = \v\ 2 V x v 
because any solution of flHOj) satisfies fl83|) . Inserting v x (V x v) x v = \v | 2 V x u into 
(IHTT) yields a cancellation and (jSTj) becomes 

V x v = ±v x (v ■ V)v. (82) 
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This is perhaps the most concise cubic form of the curl equation for v. 

The curl equation (152]) paired with the divergence equation (154"]) forms a closed 
system of first-order vector PDE for v. By construction, this set of equations for v is 
equivalent to the scalar equations ([3]), for the appropriate choice of sign. To make them 
well posed the vector equations need to be supplemented by asymptotic conditions for 
v at spatial infinity, which as before we take to be vanishing in agreement with the 
monopole asymptotics of u at spatial infinity. 

5.4-1 Spin-off: a-priori differential identities for solutions 

The cubic version of the v problem reveals two a-priori differential identities which are 
satisfied by any solution of the v problem. Namely, dotting (150]) with v yields 

v(V xv)=0 (83) 

wherever V x v is defined; and dotting (1501 with V\v\ 2 we see that 

V|f| 2 -Vxt; = (84) 

for any solution. Clearly, neither (1531 nor ([51)1 are generally true for arbitrary vector 
fields v. 

5-4-2 The Helmholtz decomposition for the cubic version 

The Helmholtz decomposition v = v g + v c leads to analogous conclusions for the cubic 
version of the v problem. We find the previously obtained linear equation (1771) for v g 
(together with V x v g = 0), while the nonlinear equation for v c (given v g ) becomes 

V x v c = ±(v g + v c ) x (v g + v c ) ■ V{v g + v c ) (85) 

(together with V • v c = 0). This closed set of first-order vector equations for the vector 
fields v g and v c is supplemented by the asymptotic conditions that v g and v c vanish at 
spatial infinity. 

As announced, given the gradient field (179|) . the remaining equation (155]) now be- 
comes a closed vector equation for the solenoidal field v c . It remains to solve equation 
(1531 with Vg given by (j79]) . 

5-4-3 Cubic version of the hierarchical v series solution for small H 

In this subsection we apply our solution method to the cubic v problem when H is 

small. We begin by deriving the cubic analog of the hierarchy of linear PDE. Since 
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H £ Cq fl L 1 is small in the sense that H(x) = eH Q (x), where H (x) is a fixed function 
and e £ R+ is chosen so small that |-H"(ar)| <C 1 uniformly in x, as before we find 
v g = 0(e), with v g given in f!79|) . Thus we write t> 9 = ei/ 1 ), and for the solenoidal 
part of v we make the power series Ansatz v c = e Pl v ^ + e P2 v^ + e P3 i/ Ps ) + . . . with 
1 < p\ < P2 < Pz < ■ ■ ■ , and with each v independent of e. Inserting this Ansatz 
into the cubic vector PDE (|85)) and identifying p\ with the smallest power on the RHS, 
Pi with the next-to-smallest power, and so on, we find recursively that p\ = 3, then 
P2 = 5, and pk = 2k + 1 for general k £ N. Moreover, for v = v g + v c with v g = ev^ 1 ' 
and v c = e 3 v^ + e 5 v^ + t 7 v^ + . . . to be a solution of the pair of equations (jSH), (|85l) . 
each t>( 2A:+1 ) for £ N has to satisfy 

V x */ 2fc+1 ) = ± vi2h+1) x ( vi2t+1) ■ V)^ (2J+1) (86) 

h+i+j= 
k-l 

supplemented by the solenoidality condition 

V • v {2k+1) = 0. (87) 

Supposing that known for all I < k — 1, then (1861) . (1871) is a pair of linear first- 

order PDE for y( 2k+1 \ with vanishing conditions at spatial infinity for v^ 2k+1 \ Now t)' 1 ' 
is known, and so, by induction, it follows that (1861) . (JHTj) successively determine u' 3 ', 
then d™, and so on — provided that each equation in this formal linear hierarchy of 
equations is solvable (in C Q ' a )\ If we ignore for a moment that we already know that our 
formal series derived from the original version of the v problem converges absolutely to a 
classical solution, from which the solvability of the cubic hierarchy of linear equations 
follows as a corollary, then the solvability of the cubic hierarchy of linear equations 
is not at all obvious but needs to be verified. Since Ylh+i+j=k-i v^ 2h+1 ^ x (ti( 2l+1 ) ■ 
Vy 2i+1) £ C Q n L 1 if ^ 2m ) £ Cl' a n L 2 for alU < k - 1, it remains to be verified 
that the collected terms on the RHS of (I8"B"1) have vanishing divergence at each order k. 
With some effort one can show term by term that (l§TI) is identical to (I6"8"|) . fIBTJl) . For 
instance, based on the identity V(a • b) = (a • V)6 + [b ■ V)a + axVx6 + 6xVxa 
for any two vector fields a and b it is readily shown that 

-v® x ( V W . vy 1 ) = iv|^ (1) | 2 x = V x (H^lV^), (88) 

proving the equality of the RHS of (157)) and that of (1551) . (15"5)) when = 1. However, 
the procedure of proving equality term by term soon gets very complicated, and an 
inductive argument is needed, instead. 
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5.4-4 Solvability of the cubic hierarchy 
Proceeding by induction, we will now show that, given v g 
,(2k+i) yfaifa vanishes at infinity and solves 



ev^ 1 ', for each fc£N there 



is a v K 



V x i/ 2fc+1 ) 



± j2 v(2h+1) x ( v(2l+1) ■ vy 2j+i) . 

h+i+j=k-l 
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Since the formulas become rather long, we switch to the shorter notation v^ for t>( 2fc+1 ) 
for all k = 0,1, 2, 3, this mildly obscures the power of e to which the terms belong, 
but shortens the length of the formulas considerably. 

Step one is easily disposed of by checking explicitly that v g x iy g -V)v g is divergence- 
free for any gradient field v g . Indeed, since V x v g = for any gradient field v g , first of 

-V\v I 2 

2 V I U 9 1 ' 



all the well-known identity (v-V)v = |V|t> | 2 — v x (Vxt>) reduces to (v g -W)v g 

so that v g x (vg-'V)Vg = — V : 



and second, v g x 



|vki 2 



-Vx 



2\ V 9\ 



2\ V 9\ 



'91 



for either sign of "±" there is a 



• V)f <-°>. Moreover, since v g is given by ( !79|) and 
V)v g e Co' a nL\ and so in particular there is a 



is a curl, i.e. divergence- free. Hence, given ev 
i>M satisfying V x v M = ±t>[°] x (t> 
H e Cq fl L 1 , we have that u s x [v. 

solution t>M of V x t>M = ±v^ x (v^ ■ V)u^ which vanishes at infinity. 

As to the induction step, suppose that for all k <m, and either sign of "±," there 
is a solution v ^ of (159]) which vanishes at spatial infinity. We now show that then 
also Y^h+i+j=m v ^ x ' V)v^ is divergence free and vanishes sufficiently rapidly at 
infinity. 

We will need two equalities. Taking the dot product of both sides of ( 1591 with 
and summing for I + k = m! with m! < m gives 



« w -v 



X V 



[k] 



k+l=m' 



=m' — k — 1 

± ^ uM • |> x (uW • V)^' ] 



M-i+j+i 
=m'-l 



(90) 
(91) 

= 0, (92) 

where the last equality follows by the anti-symmetry of the triple product a ■ {b x c); 
(1921) is the order m expansion analogue of (153]) . 

Similarly, dotting both sides of ( 159]) with (v^-V)v^ and summing over k+l+p = m 
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gives 



V x v [k] ■ ((v [l] ■ V)v [p] ) = (93) 

fc+/+n=m 



fc + Z+n h-\-i-\-j 
— m =k — 1 



^ (#1 x (uW • V)^ 1 ) • ((v® ■ V)w [pl ) = 0. (95) 

h+i+j +l+n=m— 1 



Now note that 



^ x (±V(*; [il • w b ' ] ) - x (V x uW) 



h+i+j=m 

( - V x • v W)vW) + %(v® ■ v®)V x V W 

h+i+j=m 



-(w [/tl • V x + (v [h] ■ v [j] )V x 



where we used that S h+i+i=TO (u [fc] • V x v [j] )v^ ] = 0, by QS2}. And so, 



v • v[h] x • v )^ y] = 1 v ^ [/l] ■ uM ) ■ v * '' 



. J j] 

h-\-i+j=m h+i+j=m 

[j] 



X V L 



[((«W-V)r;W + (t;M.V)t;W) -V 
+ (V' l] x V x «W + x V x i^)) • V x «W 



h+i+j=m 



By ([95D, E/»+t+i=m ((^ ■ V )^ H + ( vW ' • V x W M = 0. Moreover, by the 

anti-symmetry of the triple product a ■ (b x c) we have 

( V W x (V x W M) + uW x (V x u^)) • V x W W = o. (96) 
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Thus Ylh+i+3=m v x ( v ' is divergence free, as claimed. 

The upshot is: for all k G N the pair of equations ( 186]) . flHTj) has a unique Cq° fl L 2 
solution given by 

v W°+i)( x) = ±±J2f ( v(2h+1) x (^ (2i+1) • Vy 2i+1) )(y) x T^r^V (97) 

h+i+j J \ x y\ 

= k-l 

In addition we have ev^ = v g given by (179]) . 
5.^.5 Convergence of the cubic v series 

Multiplying ([97D by e 2fc+1 and summing over k and adding et/ 1 ) = v g given by ( 179]) 
yields a formal series solution to the cubic v problem, viz. 

« W = ±± E ^ 2,=+1 E V x / ^ X( ", _ , V) " ^dV (98) 

fc=o h+i+j= J \ y\ 

k-l 

We now have to address the convergence of this formal series solution. 

Since the combinatorial structure of the RHS of ( 19?]) is considerably simpler than 
that of ( )68|) . ( )69|) . it would seem that a convergence proof is more readily forthcoming 
than our previous proof. Curiously, the obviously simpler series expansion for the 
cubic version of the v problem has not at all yielded to our attempts of proving its 
convergence directly, i.e. without recourse to the original version of the first-order 
vector problem with its more complicated nonlinearity. So we are finally forced to 
recall the origin of the cubic version of the v problem to conclude that (198]) converges 
absolutely to a classical solution for small enough e. 

Summary: Except for the first few low-order terms, the series solution to our first-order 
v vector problem and its spin-off, the solutions to the first-order w vector problem 
and the original second-order scalar u problem, soon involves terms which look more 
and more unwieldy. The cubic version of the v problem offers relief, by having to 
evaluate fewer integrals. But there is also a caveat: it would be good to have a simple 
convergence proof directly for the cubic hierarchy. 

6 Related quasi-linear problems in divergence form 

Our power series technique of solving the prescribed mean-curvature problem for graphs 
over M n can handle more general quasilinear problems in divergence form. In particular, 
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let f(s) = 1 + as + bs 2 + • • • be given and analytic about s = 0, then the equation 



with p G Cg R L 1 small enough and asymptotic monopole condition for u, can be 
solved with the same solution techniques developed here for the scalar prescribed 
mean-curvature equation. In particular, such type of problems occur in the theory 
of stationary compressible fluid flows, see |5] ETJ [321 [33], and other versions of nonlin- 
ear electrostatics [T5] . 
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A The electrostatic analog 



In Maxwell's classical electromagnetic theory, Coulomb's law states that an electro- 
static charge density p in M 3 is the source for the electric displacement density field D, 

V • D = 4vrp Coulomb's law , (100) 

while Faraday's law says that the electric field strength E, when stationary, is curl-free: 

V x E = Faraday's law (stationary) . (101) 

The two fields E and D need to be linked by an "aether law," with the help of which 
one can eliminate either D or E and obtain a closed nonlinear set of first-order PDE 
for the remaining field. Note that "aether" is used here merely as a shorthand for 
"electromagnetic vacuum;" historically Maxwell of course thought of "the classical 
aether." 

Of particular relevance for the prescribed mean-curvature problem is the aether law 
proposed in the 1930s by Max Born j7], which reads 

E 

D = - Born's law (102) 

^\-^\E\ 2 

which can be inverted to yield 

E = — — Born's law (reverse). (103) 

y/l + P*\D\ 2 

Here, (3 G (0, oo) is a hypothetical new constant of nature (in the dimensionless notation 
of [23]). The limit (3 — > of Born's law is Maxwell's law of the "pure aether", D = E. 

The system of electrostatic Maxwell-Born (-Infeld) equations is clearly equivalent to 
the presribed mean-curvature problem in Minkowski spacetime IR 1,3 . Replacing /3 4 — > 
—(3 4 in (11021) yields a "pseudo Born's law," in which case the electrostatic Maxwell- 
Born(-Infeld) equations are equivalent to the presribed mean-curvature problem in 
Euclidean space M 4 . 

Remarks: Born's aether law fll02j) is only the electrostatic special case of his elec- 
tromagnetic aether law (replace \E\ 2 by \E\ 2 — \B\ 2 , and add a counterpart for the 
magnetic induction field B and field strength H), see and also App. VI in [9]. 

Born's electromagnetic aether law was subsequently generalized by Born and Infeld 
[TO] to a much more interesting law, yet in the electrostatic limit all these aether laws 
coincide. Also higher-dimensional generalizations of Maxwell-Born-Infeld theory with 
geometric significance exist; see, for instance, 
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